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If} . By averaging over an ensemble of field configurations, a classical field theory can 

display many of the characteristics of quantum field theory, including Lorentz invari- 
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Abstract 



CO ■ ance, a loop expansion, and renormalization effects. There is additional freedom in how 

>• ■ 

\q • the ensemble is chosen. When the field mode amplitudes have a Gaussian distribution, 

co 

CN 

CNJ ■ tween the classical and quantum theories. When the mode amplitudes are fixed, or 

have a nongaussian distribution, the quartic and higher correlations among the free 
fields are modified, seemingly in a nonlocal way. We show how this in turn affects the 
perturbative expansion. We focus on A</> 4 theory in 1 + f dimensions and use lattice 
O ■ simulations to augment our study. We give examples of how these nonlocal correla- 

tions induce behavior more similar to quantum field theory, at both weak and strong 
coupling. 
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1 Classical fields masquerading as quantum fields 

We begin by considering free field configurations of a classical field theory having the same 
energy spectrum as the vacuum fluctuations of quantum field theory. From an ensemble of 
such configurations the n-point correlation functions can be obtained, and we will discuss 
some dependence that these functions have on the choice of ensemble. The relation to 
quantum field theory will become more clear in later sections, where we discuss how the free 
configurations are to be perturbatively corrected to account for interactions. The choice of 
ensemble that we discuss here will have a nontrivial effect on this perturbative expansion. 
For a classical field theory with Lagrangian 

£ = - imV* - ~A0 4 , (1) 
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we may consider a field configuration consisting of a sum over the free modes, 

0o (#) = 2^ a P cos(co> p t + p ■ x + Op). (2) 



These are solutions for A = and lu p = y p 2 + m 2 . We will typically work in a finite volume 
V and so have a discrete set of modes. We constrain the value of the free Hamiltonian for 
these configurations to be 



and thus fix the amplitudes as 



H,. 
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Now consider an ensemble of such configurations, differing only through the choice of 
random phases P which we take to be uniformly distributed. In the limit of a sufficiently 
large number of configurations we can obtain the expectation value of a product of fields by 
simply averaging over the phases. This gives 12] 

((po(x)My)) = D °( x - y) = Yl o — v cos (p( x ~ f ))> with p = p)- ( 5 ) 
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In the V — » oo limit we could write this as 

(Mx)My)) = J ^G%( P )e-^-y\ (6) 

G%( P ) = hn5(p 2 -m% (7) 

thus showing that the \hw spectrum is Lorentz invariant in this limit. 

A more conventional choice of ensemble would have the mode amplitudes also varying 
from one configuration to the next, in such a way that only the average value of the energy 
of each mode is \hjJ v . Then the amplitudes can have a Gaussian probability distribution, 

p(a p ) = -j± exp(-^f), a\ = h/(Vu p ), (8) 

while treating the a p as nonnegative polar variables. The result (JSJ) for the 2-point function 
again follows. This choice bears a resemblance to the state vectors of a quantum field theory 
when written as wave functionals Q[4>], since the solution to the free functional Schrodinger 
equation is of course a Gaussian. 1 

The choice of ensemble affects the quartic and higher free n-point functions. They are 
obtained as before by taking the ensemble average of products of fields, which in turn are 
sums of modes. For example the free 4-point function in the Gaussian case is as expected, 

((^o(x 1 )(f)o(x2)(f)o(x 3 )(f)o(x i )) G = 

D°( Xl - x 2 )D°(x 3 - x 4 ) + D°( Xl - x 3 )D°(x 2 - x A ) + D°(x 1 - x 4 )D°(x 2 - x 3 ). (9) 



1 Gaussian distributions also arise for classical fields in thermal equilibrium, but the ensembles we are 
considering are very far from thermal. 



When the amplitudes do not vary between configurations, the fixed amplitude case, there 
are "unusual" contributions in the 4-point function due to terms where the product of modes 
involves only a single mode, so the four factors share a single phase and a single amplitude. 
The average is now over the phase only, and this changes the result to 

(M x i)M x 2)M x 3)M x 4))e = (M x i)M x 2)M x 3)M x -i))G 

- -A(x 1 ,x 2 ;x 3 ,x 4 ) - -A(x 1 ,x 3 ;x 2 ,x i ) - -A(x 1 ,x 4 ;x 2 ,x 3 ), (10) 
A(x 1 ,x 2 ;x 3 ,xa) = \ cos(p(xi - x 2 )) cos(p(x 3 -ar 4 )). (11) 

Alternatively we may write this as 

(0o(^l)0o(^2)0o(^3)0o(^4))e = 

[[D°(xi - x 2 )D°(x 3 - x 4 ) + D°(x\ - x 3 )D°(x 2 - x 4 ) + D°( Xl - x 4 )D°(x 2 - x 3 )]If 2) 

[[..]] i indicates that we must reduce the terms with a common p, after inserting (j^J for D°, 
by a factor of 1/2. 

In the case of the free 6-point function, contributions with two of the momenta the same 
are again reduced by 1/2, while the terms with all three of the momentum the same are 
reduced by a factor of 1/6. More generally, contributions to free correlation functions with 
n momenta the same are reduced by a factor of l/n\. 

Of course one could engineer a probability distribution p(a p ) to produce quite arbitrary 
changes in the coincident momentum behavior of free correlation functions. In terms of 
the second moment of the distribution, cr p , if the (2n)th moment for n > 1 is C n <Tp n , then 
the contributions with n momenta coincident at p are changed by a factor of C n /n\. The 
Gaussian and fixed amplitude cases have C n = n\ and C n = 1 respectively, and it can be 
seen that the fixed amplitude case minimizes the C n . A simple family of distributions is 

Op" 1 da 2 v 

p(ap) oc -Z-j- exp(-— — ), < d < oo (13) 

where Gaussian and fixed amplitude cases have d = 2 and d — > oo respectively. 

Thus we see that the freedom in the choice of ensemble modifies the free higher n-point 
functions, while holding the free 2-point function fixed. We notice that A in (|11|) does 
not factorize and it does not vanish for large spacelike separations of pairs of points. These 
modifications then violate cluster decomposition and thus lie outside the realm of correlations 
produced by a local quantum field theory. In this sense they could be considered nonlocal, 
even though we see that they can occur naturally in the free field configurations of a local 
classical field theory. We will refer to them as nonlocal for lack of a better description. As 
a measure of the nonlocality in the 4-point function we have 

P = 2 - C 2 = ^-1, (14) 



so that [[..]]i in (fT2|) is replaced by [[••]] 1-/3/2- 

In the classical field theory these free correlations form the basis of the perturbative 
expansion, and thus the first question is whether this expansion can be significantly affected 
by the nonlocal correlations. The perturbative expansion is a loop expansion involving 
momentum loop integrations, and it might appear that the nonlocal correlations would 
have vanishing effect in the infinite volume limit since the correction terms involve fewer 
momentum sums. On the other hand the loop integrands may display singular behavior in 
regions of momentum space where these corrections contribute. We will show, both through 
a direct sample calculation and by lattice simulation, that the nonlocal correlations do in 
fact alter the perturbative expansion in significant ways. 

The Gaussian classical theory is known to deviate from the behavior of the corresponding 
interacting quantum field theory, and we will remind the reader of the differences at the 
perturbative level. In certain situations of high temperature and/or high occupation number, 
these differences may be controlled and the classical theory used as an approximation scheme 
But at zero temperature the departures from the quantum theory that are evident for 
the Gaussian classical theory could be taken to define what truly quantum phenomena are. 
Thus it is of interest to compare these truly nonclassical effects of the quantum theory 
to the nonlocal effects of the nongaussian classical theory. We will find that the nonlocal 
correlations of the free fields can act to bring the perturbative expansion of the classical 
theory closer to that of the quantum theory. This helps to explain the findings of E] • 

In the next section we will develop the Gaussian classical theory, highlighting both its 
resemblance to quantum field theory and its qualitative differences. Section 3 will pinpoint 
the perturbative difference, thus making clear what is missing relative to the quantum theory. 
In section 4 we consider a 2-loop example where we identify the missing piece, and we compare 
it to a similar contribution coming from the nonlocal correlations in the nongaussian classical 
theory. In sections 5 and 6 we use lattice simulations of the classical theory to investigate this 
more fully, and we find further evidence of quantum-like perturbative effects. In particular 
damping and thermalization effects are greatly reduced. The simulations can also be used 
at strong coupling, and in section 7 we will see quantum-like critical behavior emerging in 
the nongaussian classical theory. We end with some comments and open questions in the 
final section. 

2 The Gaussian classical theory 

It is possible to develop a perturbative expansion by writing the field equation as an integral 
equation and then iterating |2j. Instead we shall make use of a path integral that does 
nothing but enforce the classical field equations, and then obtain the perturbative expansion 
from the path integral in the usual way. The integration over the fields in the path integral 
will produce the sum over the ensemble of classical field configurations. 



The expectation values of products of fields in some ensemble of configurations can be 
expressed as follows, 

(0(xi)0(x 2 )...) = / ^xp[^)]^ii)^(i 2 )...e^^ 4 ^ fel (15) 



where tj < {t\,t2, ■■■} < tf. (A similar construction appears for example in j^J.) The 
Lagrangian C x is linear in the auxiliary field x 

C x = d^ X - m 2 </> X - A0 3 x, (16) 

and we may require that X ( x ) vanishes at times i» and tf. The functional integral over x 
gives the constraint 

Y[ 2ttM [(□ + m 2 )0(r, x) + A0 3 (r, x)] (17) 

T,X 

for tj < r < tf. In this way the field <fi is constrained to satisfy the interacting classical 
field equation. Notice that a h~ l has been inserted in the action, even though it has no 
significance at this stage. 

p[4>(xi)], where X{ is the coordinate on the initial time-slice t = ti, selects the configu- 
rations in our ensemble. Since we are interested in a perturbation theory it is sufficient to 
describe this ensemble in terms of the free field modes, and in this section we are making 
the Gaussian choice. These configurations are acting as the initial conditions for evolution 
by the full field equations. 

Now the point of this construction is to leave the x field in the theory to develop the 
perturbation theory. The Lagrangian specifies the <ftx and x0 propagators up to boundary 
conditions, and since we are studying a classical theory it is appropriate to choose the 
retarded and advanced Greens functions, 

r<o ^ r o ^ / 18 \ 

<t>X ~ ~~2 2 i • X4> ~ ~~2 2 : " V 

With this choice we see diagrammatically that the x en d °f a li ne always appears at a later 
time than the end. But each vertex has only one x an d three 0's, so closed loops of G^ x 
or G X( p lines do not form. The resulting tree graph expansion is the expected result for a 
classical theory. 

This would be all in the absence of background fields, but due to our ensemble of field 
configurations we also have G^Jp) — hn5(p 2 — m 2 ). The h appearing here, originating in 
the mode amplitudes, does have significance. Now the branches of the various trees can be 
interconnected with G^ lines, and a nontrivial loop expansion emerges. In fact the same 
relation between the number of loops in a diagram and the power of h occurs in this classical 
theory as it does in quantum field theory. The fr 1 in the action was not necessary for this 
but it makes this result easier to see, since then factors of h are associated with propagators 
and vertices in the usual way. 



The free 2-point functions are also such that = and 

XX 

G% = \e{ Po )(Gl x -G^). (19) 

We can make some statements about the full 2-point functions (denoted by G) and the 2PI 
self-energies (denoted by II). Diagramatically it can be checked that = 0, and similarly 
G xx = continues to hold in the interacting theory. Then the full G$ x and G x $ satisfy 

G<t>x = Gd> x + ^iJl^Gfa (20) 



G X 4> = G X 4> + G^UfaGxtf, (21) 



and thus 



G * = jfi-mi-lC+iilU + epoy , '^ ElI '" in< (22) 

ih 

G x<t> = — 2 2 tt TFi—, V' im 4>x = ^r + HU (23) 

p — 77i — II r — i(Ui + epo) 

These G's have the same analytic structure as in the free theory, so IT is an odd function of 
p with Ilj(po) > for p > 0. 

Of more interest is the full 2-point function G^ for the original scalar field. The sum- 
mation of the diagrams involves various classes of diagrams. 

= G^ + G^n^G^ + GfaH-x^G^ + G^U^G^Il^G^ + G^ x U xx G xll> (24) 
= G<^ x n xx G x ^ (25) 

The cancellations among the various classes follow from the expressions for GjL,, G^ and 
G x< f,. We can compare this to the result of using ([TO*]) as a relation among the full 2-point 
functions, which gives 

^ " (p2_ m 2_ IIr )2 + II 2 ^ 

= G^ipoWiG^ (27) 

Comparison of (|2*5|l and (|2*7|) would imply the relation II XX = e(p )IIj. This is in fact true, 
as can be checked diagramatically. 

It turns out that IT (or IT XX ) is nonvanishing on mass shell, in which case we have the 
finite width (and shifted mass due to Il r ) form of a pseudoparticle correlator. The emergence 
of a pseudoparticle is typical in classical field theory applications, in thermal systems for 
example, and the pseudoparticle width is often referred to as plasmon or Landau damping. 
A sampling of references to this and related applications of classical field theories is given in 

m 

This pseudoparticle distinguishes the classical theory from the particle description of a 
quantum field theory. Nevertheless there is a loop expansion characterized by powers of h, 



and a renormalization of the parameters m and A in C x similar to quantum field theory. In 
addition to the XX amplitude, the classical theory will perturbatively generate amputated 
amplitudes of the form <fr r x s with r + s even. We have seen that a <f)<j) self-energy is not 
generated and more generally amplitudes are not generated unless s > 1. A real contribution 
to the (f>x 3 amplitude is expected, and in the next section we will find this amplitude to be 
of special interest. 

[Another implication is that C x is not corrected by a constant term, a cosmological 
constant. Here it is useful again to emphasize the difference between the Lagrangian C x 
appearing in the path integral and the original classical Lagrangian C in It is the latter 
the defines the energy momentum tensor, whose the expectation value can be calculated 
perturbatively using the path integral involving C x . A contribution to the cosmological 
constant could be obtained in this way] 



3 A small change and the true quantum theory 

Now let us consider a new theory by adding a <ftx 3 interaction to our previous theory defined 
by C x in (jl5H16|) . where we continue to assume the Gaussian ensemble. Denoting the new 
Lagrangian by C a , 

C a = C x + a^ X 3 , (28) 

it is clear that the direct link to the classical evolution via integrating out the x field has 
been lost. In fact this additional interaction term is all that separates the classical theory 
from the full quantum field theory. (This observation [3J of the relation between classical and 
quantum field theories originates in studies of classical fields used to model quantum fields 
at high temperature or high occupation number.) This may be somewhat surprising, given 
that we have just observed that this amplitude is generated at some level in the classical 
theory. 

The following redefinition of fields 

X= J -(</>+-<!>-), = ^(0+ + 0-), (29) 

V OL I 

yields 

£ a = y^ [£(</>+)-£(</>_)], (30) 

where £ is the original classical Lagrangian in Thus C is now appearing in the path 
integral itself. If a = 1 the overall factor in the action is the usual one for a quantum field 
theory, given the h~ l factor already in the action of ()15|). 

In this case of a = 1 the transformed versions of the original free propagators take the 



following form. 



G l + = iAr ( 31 ) 

= 2 ^ . (32) 

G\ = 2h7r6(p )5(p 2 - m 2 ) (33) 
G°_ = 2h7r6{-p )5{p 2 - m 2 ) (34) 

G+ + is the Feynman propagator. We have now arrived at the Schwinger-Keldysh formalism 
for the calculation of "in-in" matrix elements in quantum field theory jHj. This formalism 
is complementary to the more standard formalism of quantum field theory, concerned with 
the calculation of "in-out" matrix elements. 

For example jS] consider the full G ++ 2-point function 

G ++ (x u x 2 ) oc [ dtf>+dtf>. p[0 + (x l )]0 + (x 1 )0 + (x 2 )e*4 /d4 ^ + )-£W-)) (35) 



The original constraints on the x fi e ld now read 0_(xj) = <fi + (xi) and 0_(x/-) = <p + (xf). 
Other than this the <p + and 0_ path integrals have decoupled, and each is a path integral 
that can be identified with a "in-out" transition amplitude or correlation function of quantum 
field theory. Thus 

G ++ = y^(in, vac|out, s)(out, s|T0(xi)0(x 2 )|in, vac) (36) 

s 

= (in, vac|T0(xi)0(x 2 )|in, vac) (37) 

G , G + - and G_ + have analogous expressions. The sum over "out" states corresponds 

to the fact that no constraint has been placed on <f){xf) in the path integral. And p[(p(xi)} 
has given us an explicit representation of what we mean by the "in" vacuum state; for the 
Gaussian choice we have completed the equivalence to quantum field theory. 

Thus the modification of the theory through the addition of the <px 3 term has replaced 
the classical correlation functions by their quantum counterparts. We see that the whole 
essence of quantum field theory (for a scalar field) lies in a particular value for the coefficient 
of the <frx 3 term. The essence is not the h, or the loop expansion, or renormalization, since 
all this appears in the classical theory. And as the coefficient a varies from to 1, the theory 
extrapolates continuously from classical to quantum. 



4 Quantum behavior at 2-loops 

In this section we first isolate a truly quantum perturbative contribution to the self-energy, 
a contribution that is due to the <px 3 term. We will see explicitly that this contribution for 
a = 1 has the effect of cancelling the classical contribution to the imaginary self-energy IT on 



classical 



quantum (oc a) 



U 4>X 




^XX 



Figure 1: The 2- loop self-energy diagrams are split into classical and quantum pieces. The 
solid lines are G° x or and the dashed lines are G^. 

mass shell. Then the width of the pseudoparticle goes to zero, and we end up with a stable 
particle as described by quantum field theory. But the main point here will be to continue 
our discussion of section 1, and obtain the nonlocal corrections to the same diagram that 
can occur in the classical theory. 

A truly quantum contribution to the self-energy first appears at the 2-loop level. We 
return to the (0, x) basis where a 4>x 3 term generates loops of G° x or G^ lines, thus allowing 
more than one or G^ line to connect two vertices. We show the resulting classical and 
quantum pieces of the sunset diagram in Fig. (JJJ). We have seen that the IT XX diagrams are 
just the imaginary parts of the H^ x ones, so we need only consider the top two diagrams. 

We restrict ourselves to 1 + 1 dimensions as we do in the lattice simulations to follow, 
and consider the classical 2-loop contribution to H^, the first diagram in Fig. (JTJ). Without 
Wick rotating and in finite volume V = L this is 

n 2-ioo P (k) = 6! 2 A 2 [ dp dq tx5 (pg - p\ - 7ti 2 )tt5 (gg - q\ - m 2 ) 
classical ) 2L 2 ^ J 2tt 2tt (fc - Po ~ Qo + ~ (fci ~ Pi ~ qi) 2 ~ m 2 ' 



Pi,<n 



After the temporal momentum integrals are done we are left with sums over the spatial 
momenta p\ and q\ whose discrete values are multiples of 2ir/L. These sums are infrared 
dominated. They are evaluated numerically while keeping the e in the retarded propagator 
small but nonvanishing. Of interest is the on-shell self-energy, obtained by setting (fc , ki) = 
(m, 0). In a similar way we may consider the quantum contribution proportional to a, the 
second diagram in Fig. (JJJ. 

Combining the results for the real and imaginary parts we have 

U 2 r - loop {m) oc l-ct/3, (39) 
nf oop (m) oc I -a. (40) 

Thus in the quantum theory defined by a = 1 we find that n^" loop (m) is 2/3 of its size in 
the classical theory. And n 2_loop (m) vanishes as expected. The classical and quantum pieces 



of n^" loop (m) are each finite and quite insensitive to how the V — > oo and e — > limits are 
taken. For U 2 ~ loop (m) on the other hand there is a weak divergence in each piece as e — > 0, 
only canceling when a — 1. We note that this divergence is an artifact of the perturbative 
calculation, since the full IT, (m) is certainly finite in the classical theory. We will see this in 
the simulations. Although we will not pursue it here, one could consider replacing the free 
propagators by the full propagators and G^ x in the 2-loop calculation. 

Now that we have isolated the "truly quantum" from the classical, we may consider the 
effect of making another choice for the ensemble in the classical theory. We return to the 
classical theory defined by a = 0, but now account for a modified perturbative expansion 
corresponding to the different choice of p[4>{xj)\ appearing in (j!5|) . According to section 1, 
as we depart from the Gaussian prescription the free n-point functions are no longer simply 
expressed in terms of the free 2-point function. We parametrized the nonlocality in the 4- 
point function by (3, where (3 = and (3 = 1 correspond to the Gaussian and fixed amplitude 
cases respectively. The question is how a nonvanishing (3 affects the 2-loop calculation in 

(EH). 

The Feynman rules for the two CrjL lines in the diagram arose by writing the D° in (JHJ) 

as 

d°(x -y) = Y,vf tr G > )e ~" ( ^ (41) 
pi 

We notice that relative to (jSJ) this incorporates some inconsequential changes of sign of p\ 
in the terms that have negative p$. But to identify which contributions in ()38|) need to be 
modified, we should undo these changes of sign by replacing pi inside the sum in (|41j) by 
sign(po)pi, and the same for p\ and q± in Then the rule from section 1 is to reduce 

the terms with pi = qi in the discrete momentum sums by a factor of (1 — (3/2). The net 
result is easier to describe in terms of the two 2- vectors p^ and q^, now with the sign factors 
included in their definition. The rule is to reduce by (1 — (3/2) the contributions where the 
2- momenta of the two G^ lines satisfy either p^ = g M or p^ = —q^. 2 

But this is not all, since what is directly measurable is the 2-point function rather than 
the self-energy. And one of the two external lines of the sunset diagram contribution to the 
2-point function must be a G^ line (second diagram in Fig. (J2J). The same is true at 1-loop 
(first diagram in Fig. Q), which in this case implies that there are two G^ lines in that 
diagram. Thus there is also a nonlocal correction to the 1-loop contribution to the 2-point 
function. For an external line at vanishing spatial momentum, the correction involves the 
zero momentum contribution to the loop. This then is a factor of (1 — (3/2) smaller than 
usual while all the other modes contribute in the loop as usual. It can easily be checked that 
this nonlocal correction at 1-loop will vanish in the V — > oo limit for fixed mass. 

Returning to the sunset diagram, the external G^ line means that there are three G^ 
lines in this contribution to the 2-point function. Thus when either of the internal G^ lines 

2 Wc notice the Lorentz invariancc of this rule. 



Figure 2: 1-loop and 2-loop contributions to the 2-point function. 



has vanishing spatial momentum, these contributions should also be reduced by (1 — (3/2). 
Unlike the 1-loop case, we find that these corrections, along with the previous p^ = or 
p^ = —q^ corrections, do not vanish in the V — > oo limit. 3 They produce a finite correction to 
n^" loop (m) as long as the limit e — > is taken before the limit V — > oo (in the opposite order 
it vanishes). More precisely e must be sufficiently small compared to lattice spacing 27r/L 
of the momentum space lattice used to evaluate the sums. The finite correction illustrates 
how the singularity structure of momentum loop integrals, arising from e — * 0, can offset the 
naive volume suppression of the correction. 
The actual result from numerical analysis is 

U 2 r Aoop (m) oc 1-/3/2, (42) 
Uf loop (m) oc 1 -(3/2. (43) 

Each of the three ways of having pairwise coincident momentum in the three lines 
contribute equally in il^ loop (m), whereas only the case of coincident momentum on the two 
internal lines contributes in n?~ loop (m). But since each result is proportional to 1 — (3/2, 
this shows that only the coincident momentum terms contribute; the contributions without 
coincident momenta vanish! 

We see that n,^" loop (m) can be reduced from the Gaussian classical value as in the quantum 
theory. ITj~ loop (m) is still plagued with the overall e — > divergence, so result (|4H|) is not 
particularly useful. To see how this quantity is actually behaving we turn to the lattice 
simulations of the classical theory (a = 0). We will be comparing the Gaussian (3 = and 
the fixed amplitude (3 = 1 cases. 

5 Lattice simulations 

We wish to obtain both the real and imaginary parts of the self-energy through a direct 
lattice study of interacting classical fields in 1 + 1 dimensions. The details of the simulations 
are found in E] . They are performed by numerically evolving each configuration according 
to the full field equations in real time. For lattice spacing a, the spatial size of the lattice 
is Na where N = 256 unless otherwise specified. The definition of the initial configurations 
requires a mass parameter, and this is chosen to match the physical mass as determined by the 
3 The contribution when all three spatial momenta vanish does vanish. 



simulation. The expectation values of products of fields are directly obtained in coordinate 
space. An accurate determination of the mass corrections can be more easily determined 
from time-like rather than space-like correlators, since the latter falls exponentially. We will 
consider two time-like correlators, where the first one involves the zero mode only, 

D zm (t) = J ^G^(p)\ pi=0 e-^\ (44) 

D *® = zE/t^o°) e ~^ ( 45 ) 
pi 

In the free case with G^ — > G^ = hn5{p 2 — m 2 ) we have 

D° zm (t) = ±- cos(mt), D° t (t) = \Y,^ cos ^- ( 46 ) 

pi 

When we instead use G^ from (J26|) . and use a narrow width approximation where the 
self-energy functions are replaced by their on-shell values, then D zm (t) takes the form of an 
oscillation with decaying amplitude. This is the standard manifestation of plasmon damping 
[Zj, where the plasmon damping rate is 7 = IT(m)/2m. 

We thus determine the plasmon decay rate 7 by fitting the resulting decaying form of 
D zm {t) to what emerges from the simulation. We compare 7 for the Gaussian and fixed 
amplitude cases for a range of ma. For (ma)" 1 < 5 we find that 7's are comparable. For 
larger (ma) -1 the two 7's rapidly depart from each other, and for (ma) -1 = (15,23,30) we 
find that the 7 in the fixed amplitude case is ~ (6, 28, 45) times smaller than for the Gaussian 
case. 4 This is for small coupling where the leading 2-loop contribution should dominate. 5 
Thus we once again see that the choice of ensemble has a significant affect at the 2-loop level. 
Given that 7 should vanish in the quantum limit, we see that the fixed amplitude ensemble 
is causing the lattice simulation to move towards quantum-like behavior. 

We should stress that "fixed amplitudes" refers only to how the amplitudes are chosen 
in the initial configurations; after the initial time the amplitude of each mode is free to 
evolve as the dynamics dictates. The large reduction in the plasmon decay rate indicates 
that these configurations are more stable; thermalization processes are much slower than for 
Gaussian classical simulations. When we next study the real part of the self-energy we are 
forced to use the fixed amplitude ensemble, since it is the slow thermalization that makes 
the investigation of this quantity possible. 

n^" loop (m) is more difficult to extract than n 2_loop (m), since there is also the real 1-loop 
contribution. The latter is not infrared dominated (it is log divergent in the continuum) and 
thus completely dominates the 2-loop contribution for small coupling. To isolate the 2-loop 
4 We also checked that analogous behavior occurs on a larger lattice, N = 512. 

5 These ratios are obtained by considering different couplings in the two cases that give comparable 7's, 
and then accounting for the different couplings by multiplying the ratio of these 7's by the inverse ratio of 
squared couplings. 



contribution, values of H r (m) may be extracted for a range of the effective dimensionless 
coupling A/m 2 so that both the 1- and 2- loop contributions can be fit simultaneously to 
the measured linear plus quadratic dependence. This procedure was used in jl], where its 
relation to the 2-loop gap equation was described. 

We shall adopt a complementary procedure here and will vary H rather than A/m 2 , since 
the 1-loop and 2-loop diagrams are proportional to h and ft 2 respectively. In the classical 
theory h originates in the normalization of the mode amplitudes, so we can just vary this 
overall normalization. It is a good check on the simulation to find that these two methods 
give similar results. As we vary h we also vary the bare mass mo in the field equation so as 
to hold the physical mass m fairly constant. It is the difference between m and m that is 
fit to terms linear and quadratic in h, and thus obtaining the 1- and 2-loop contributions to 
II r (m). 

For this analysis we use the 2-point function D t (t) rather than D zm (t). In 1 + 1 dimensions 
the form of D t {t) reflects the contributions from both high and low momentum; it has both 
a rapidly oscillating component and an oscillation with 1/m time scale. Since D t (t) receives 
contributions from all the modes, it is expected to be a more reliable quantity than D zm (t). 

We first consider the 1-loop mass correction, since the extracted value of this quantity 
displays a small deviation from the naively expected result, as was noted in Indeed 
this discrepancy arises from the nonlocal correlations, as we have already described in the 
previous section as a finite volume effect. The effect increases for smaller ma, and amounts at 
most to about a 3% reduction in the 1-loop diagram. The observed and calculated values of 
this discrepancy agree, and this is a further indication that the simulation sees the nonlocal 
correlations. Returning to our fit of the linear and quadratic behavior, we now input our 
calculated linear behavior and use the data to extract the quadratic behavior only. 6 

We compare the simulations to the explicit calculations of classical (Gaussian with a = 0) 
and quantum (a = 1) values of II 2 " loop /m 2 , which at A/m 2 = 0.4 are 0.022 and 0.015 
respectively. In Fig. (J3J) we display the extracted values of II 2 " loop divided by the classical 
value as a function of (ma) -1 . We see that as ma decreases that II 2 " loop gradually drops 
below the quantum value and becomes consistent with the calculated value of 1/2 from the 
previous section. 

In Fig. (jlj) we display the associated values of n 2_loop /m 2 , which was fit simultaneously to 
the same simulation data. In the context of the plasmon decay rate we have already noted 
how small this quantity is compared to its Gaussian classical value. In the latter case this 
quantity increases with increasing (ma)" 1 , but with the fixed amplitude prescription it falls 
towards the vanishing quantum value. 

6 The fact that this was not done in 0] accounts for the slightly different results there. 



0.6 



0.4 



0.2 



T 

self-ener£ 

: 1 1 | 

Mil 


;y, re£ 


a F 


art 
















(ma) 1 









5~ 10 15 20 25 30 35 



Figure 3: The real part of the 2-loop self-energy, U 2 ~ loop (m), extracted from the simulation 
and divided by the Gaussian classical result. The quantum result at 2/3 and the fixed 
amplitude result at 1/2 are also indicated. We set X/m 2 = 0.4 and let h be .25, .5, .75 and 
1. We evolve 5000 configurations for each value of h and then perform a fit. We repeat this 
10 times to obtain an estimate of the error, for each value of ma. 
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Figure 4: The imaginary part, U 2 ~ loop (m) / m 2 , extracted from the simulation. Unlike the 
previous figure, only h = 1 is used. 



6 Interpretation 



We have used lattice simulations at weak coupling to extract 2-loop effects, and have con- 
firmed that the choice of ensemble can have a significant effect, by comparing the Gaussian 
and fixed amplitude prescriptions. We then focused on the latter and found interesting and 
almost quantum-like behavior in the real and imaginary 2-loop corrections for a range of 
ma. The question now is what is the meaning of this range of ma. 

Some insight is provided by our direct evaluation of the 2-loop mass shift Il^" loop (m) in 
section 0J For the latter we found that we had to take the limit e — > before the V — > oo 
limit. If V is measured with respect to 1/m then for the simulation with a fixed number of 
lattice points, increasing V corresponds to increasing m. And e corresponds loosely to the 
measured plasmon decay rate 7. Thus if the same order of limits is to apply in the simulation 
then there is an upper bound on V or m, since 7, although small, is not exactly zero in the 
simulation. Thus we see why quantum-like behavior requires smaller m. 

What is nontrivial is the fact that 7 is as small as it is, allowing nearly quantum-like 
behavior for a range of ma. Our direct evaluation in section 0] did not anticipate this, 
suffering as it did with a e — * divergence. The real part of the self-energy was properly 
anticipated, but this basically assumed that 7 would turn out to be sufficiently small. 

We note that the interesting behavior is occurring for values of (ma) -1 that are smaller 
than where conventional finite volume effects are expected, which is when (ma) -1 becomes 
closer to N. In fact finite volume effects do appear to show up in the simulation for (ma) -1 > 
35. There the mass extracted after long evolution times shows a tendency to drift downward. 
A drift down would result in spuriously large values of Il^" loop (m) , and thus these finite volume 
effects behave oppositely to the observed departure of n 2 ~ loop (m) from the Gaussian classical 
value. 

To more directly test for conventional finite volume and discretization effects we per- 
formed a direct lattice calculation of the 2-loop graph from quantum field theory, after Wick 
rotation jlj. The loop integrations become sums over the discrete Euclidean momenta on a 
N x N lattice. 

N/2-1 

n 2ioo P iatt (fc) = _ J- G(p 1 + k,p 2 )G(q 1 ,q 2 )G(p 1 + qi ,p 2 + q 2 ), (47) 

{Pi,P2,qi,Q2}=~N/2 

G(Pl ' P2) = 4 sin 2 (7r Pl /A0 + 4 sin 2 (7rp 2 /iV) + m 2 (48) 

This yields the self-energy for space-like or vanishing momentum. Of interest is the m 
dependence of the result, and we find for N = 256 that n 2 ; loop ~ latt (0)/m 2 varies by less than 
3% for 10 < (ma) -1 < 48. We may also compare this lattice result to a non- Wick-rotated 
calculation of the same diagram in the continuum limit, along the lines of section 0] but 
evaluated at zero momentum. We find agreement to within a few percent in the same range 
of ma. Thus conventional discretization or finite volume effects are significantly smaller than 



the effects induced by the choice of ensemble. 7 

We mention a few more details of the simulation. D t (t) is measured over a time interval 
equal to half the spatial size of the lattice At = aN/2. We evolved each configuration forward 
for a total time equal to (n& + l)At, and then measured D t (t) on each of the At time intervals 
except the first. 8 ' 9 toa = 1 for the largest value of m and we scaled n& as 1/m for smaller 
values of m. By averaging D t {t) over this many time intervals for each configuration we 
reduced the noise in the extracted masses. By using smaller ua for larger m we avoided 
the phenomena observed in |2j, namely a slight upward drift in the masses observed for 
increasing evolution time (if the evolution time was increased much beyond (riA + I) At). 
The drift is more pronounced for larger mass. In fact we suspect that the enhanced drift for 
larger mass is related to the observation that the system behaves more classically for larger 
mass; greater plasmon damping implies enhanced thermalization processes, meaning that 
our field configurations are less stable for larger mass. 

7 Strong Coupling 

In this section we shall describe some evidence of quantum-like behavior occurring beyond 
the 2-loop level. At sufficiently strong coupling the quantum theory is known to have a 
broken phase where (</>) 7^ 0, with a critical line in the tuq-X plane separating the two phases. 
This is usually defined as a line of fixed g = A/mg ap where ^g ap is the 1-loop gap mass. We 
may study this phenomena in the classical theory, where a priori there is no reason for the 
onset of ((p) 7^ to occur at a coupling g similar to the quantum critical coupling. 

A direct comparison of the quantum and classical lattice theories has the complication 
that the lattice regularizations are usually different. A square Euclidean lattice is often used 
for quantum Monte Carlo simulations, whereas the classical simulations are Lorentzian with 
the discretization in the time direction finer grained than the space direction. As described 
in [1] there is a special line of constant m^ ap in the m -A plane that is the same for the 
two regularizations; thus by comparing the location of the critical coupling on this line we 
may reduce the dependence on the lattice regularization. The result from the quantum 
lattice theory for the critical coupling is g « 10 jH] , which on the special line corresponds to 
A « 50 /N 2 a 2 . 

In the classical simulation at strong coupling we no longer have the luxury of long sim- 
ulation times, but there is still a window where quantum behavior can set in before ther- 
malization effects dominate. We determine the value of (<p) by averaging over space and 

7 Discrctization could be affecting the results for (ma) -1 < 10, but this is not the main region of interest 
in this work. 

8 The first time interval allows some time for the free initial configurations to evolve. 

9 In our use of D zm {t) to determine the plasmon decay rate, discussed above, we measured D zm (t) once 
over the whole time interval n/^/S.t. Unlike D zm (t), D t (t) drastically changes its behavior for time intervals 
larger than fs 2At. 
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Figure 5: Histograms showing the relative probabilities of finding different values of ((f>) in 
the simulation, as explained in the text. The couplings are g = (7, 9, 10) from bottom to top, 
the bin size is 0.05, and 10000 configurations are used for each g. Overall normalizations 
have been adjusted to equate the values in the first bin. 

over times between Na/2 and 3Na/2, and then averaging over configurations. The mass 
parameter in the initial configuration is chosen to be the value that defines the special line, 
m gap ~ 5/iV 2 a 2 . A constant is also added to the initial configuration to provide a nonzero 
{(f>). It is adjusted to match the most probable value that is measured at the later times, 
since this value should be close to the minimum of the associated effective potential. To 
find the most probable value at later times we find the peak of the histogram obtained by 
binning the values of the field. After iterating so that the input value of ((f)) matches this 
most probable value, we then determine ((f)) as the actual average of the field at the later 
times. In the absence of a nonzero most probable value, ((f)) is taken to vanish. 

We start with the fixed amplitude prescription. In Fig. (jHJ we display three examples 
of histograms, with g = (7,9, 10), with a vanishing ((f)), a slightly nonvanishing ((f>), and a 
larger (<f)) respectively. From this we can see that the critical coupling is a little smaller than 
9. We can test for the effect of lattice size by repeating this procedure for both a smaller 
N = 128 and a larger N = 512 lattice. Here we determine the couplings that give the same 
slightly nonvanishing ((f)) as the g = 9, N = 256 case; for N = (128, 512) the corresponding 
couplings are g ~ (6.5, 10). Thus the finite size effects are diminishing for the larger lattices. 
We see that these results are quite consistent with a critical coupling of g ~ 10 as found 
in the quantum theory. (A less detailed analysis of this effect was described in j2] . In that 
reference another method to locate the critical line, involving the overlap of the time and 
space correlators, was perhaps less precise but gave similar results.) 



Now we can consider the effect of using the Gaussian distribution prescription instead, 
keeping everything else the same. Here we find that the critical behavior shifts towards 
weaker coupling; the coupling that gives the same slightly nonvanishing (<p) is now g m 3 or 
4. The clearly discernible difference between these two cases again shows the effect of the 
choice of ensemble, this time at strong coupling. And once again it is the fixed amplitude 
prescription that produces quantum-like behavior. 

8 Comments 

1) In section 1 we noted that the amount of nonlocal correlations in ensembles of classical 
fields was controlled by the probability distribution for the amplitudes. For the 4-point 
function we parametrized the departure from locality by /3, where a Gaussian distribution 
has (3 = 0. In section 3 we described a class of Gaussian theories having a new interaction 
term with coefficient a, for which a = 1 corresponded to quantum field theory. Thus the 
first class of theories is nonlocal and the second is nonclassical, and they could be thought 
of as different ways to deform the a = (3 = theory which is both local and classical. All 
these theories have h and a loop expansion. In fact the two classes of theories belong to 
a larger class of theories in which both nonlocal and nonclassical features are introduced 
in varying amounts (both a and (3 nonvanishing). But in this work we have focused on 
nonlocal effects in interacting classical theories and have shown that they can resemble the 
nonclassical effects of the corresponding local quantum field theories. These two classes of 
theories may have more in common than commonly thought. 

2) We observe that the existence of spatial dimensions are crucial for producing quantum- 
like behavior. In + 1 dimensions then there is just one site and one mode. Averaging over 
the phase and holding the amplitude fixed in this case gives results that are in no way 
quantum-like. In fact the 1-loop correction to the mass is 1/2 its quantum- mechanical value, 
as can be deduced from our discussion of the correction to the 1-loop contribution in section 
13] This can also be found through a direct perturbative solution of the + 1 dimensional 
field equation. Thus as the number of lattice points in the spatial direction becomes smaller, 
we would expect our results to move away from quantum-like behavior, and approach the 
non-quantum-mechanical behavior of the + 1 dimension limit. But in the other limit of 
fine-grained spatial dimension(s), a nonrelativistic limit of the theory could presumably be 
taken as done in quantum field theory, and in this way making an approach to something 
that has some resemblance to quantum mechanics (HJ. 

3) The stability of the classical system and quantum-like behavior are both related to 
the vanishing of the imaginary self-energy on mass shell. We found some tendency for this 
vanishing in the lattice simulation even though there was a e — > divergence in this quantity 
in perturbation theory. Can the perturbative description be improved, possibly in some 
other classical theory? 



4) The 2-loop sunset diagram turned out to be maximally sensitive to the nonlocal 
correlations of classical fields. How general is that result? 

5) In the examples we have given, the quantum-like behavior has emerged in quantities 
that are infrared dominated. Is it essential that the theory itself be infrared dominated, and 
do asymptotically free gauge theories qualify? 
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